Sharp k-space features in the order parameter within the 
interlayer pair-tunnehng mechanism of high- Tc superconductivity 
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We study the k-dependence of tlie gap function of a bilayer superconductor, using standard mean- 
field techniques applied to a two-dimensional (2D) extended Hubbard model, in the presence of 
coherent interlayer pair-tunneling and quenched coherent single-particle tunneling. The intralayer 
pairing potential thus defined is expandable in a finite number (5) of basis functions for the irre- 
ducible representations of the point-group of the perfectly square lattice, Cav This gives rise to 
a competition between s- and d-wave symmetry, as the chemical potential is increased from the 
bottom to the top of a reafistic band for most cuprates. It allows for mixed-symmetry paired state 
at temperatures below Tc, but never at Tc on a square lattice. Inclusion of the interlayer pair- 
tunneling into the effective pairing potential leads to highly non-trivial k-space structures, such as 
pronounced maxima along the Fermi line not seen in the absence of interlayer pair-tunneling. We 
show how such a gap structure evolves with temperature and with band filling, and how it affects 
various observables. In particular, a nonuniversal value of the normalized jump in the specific heat 
at Tc will be evidenced, at variance with the conventional universal BCS result. 
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I. INTRODUCTION 

The identification of the character of the asymptotic 
low-energy excitations of the the high- Tc superconduc- 
tors (HTCS) continues to present a formidable challenge 
to theorists and experimentalists in condensed matter 
physics. These excitations are presumably a key fea^ 
ture in understanding the basics of the phenomenon.El 
Although the superconducting state of the cuprates to 
a large extent in the recent past has been considered 
conventional, it is becoming increasingly clear that such 
a stajtement requires certain modifications, to say the 
least The latter statement is supported by recent ex- 
perimental findingsE'Q 

The controversy over the symmetry of the paired state 
(s- and extended s-wave vs higher order waves, partic- 
ularly c?-wave) and the coupling strength can nowadays 
be restated in more precise terms, due to the availability 
of samples with adequately pure composition and struc- 
ture, and of improvements of experimental techniques. A 
central tool in this context is. angle-resolved photoemis- 
sion spectroscopy (ARPES),eI with which one is able to 
extract, if not the phase of the superconducting order- 
parameter (OP), then at least the k-dependence of its 
modulus at various temperatures and chemical compo- 
sitions. Here, k is a wave- vector ranging over the first 
Brillouin zone (IBZ) of the appropriate inverse lattice 
for the cuprate compound under consideration. In par- 
ticular, there is a growing consensus on the occurrence of 
nodes of the OP along the — ky direction in the IBZ 



for optimally doped BiaSrajCaCuaOg (Bi2212).i However, 
some contradictory claims 13 for different samples seem to 
support, in a parallel way, the idea that the detailed k- 
space shape of the OP could be a material specific prop- 
ertiL although the location of the nodal lines may not 
be.El This is in agreement with the fact that both the 
critical temperature Tc itself and the maximum gap at 
T = change considerably from one material to another, 
as well as within a given material class, on varying the 
doping level. 

We will in this paper try to bring out a few peculiar 
details of some properties of the superconducting state 
within the interlayer pair-tunneling mechanism (ILPT), 
which seems to be almost unique to this pairing mecha- 
nism. It is at any rate becoming clear that the determina- 
tion of the location of nodal lines in the superconducting 
OP, i.e. its transformation properties under the sym- 
metry operations of the underlying lattice, by no means 
suffices to unambiguously determine the unconventional 
pairing mechanism. In this sense, the symmetry of the 
OP is perhaps not a central issue, although it certainly 
has been the focus of much research during the last few 
years. Moreover, the controversy over the symmetry of 
the OP has initiated some of the most sophistipiai^ ex- 
periments in condensed matter physics to date.u~li3 

In this paper, we shall mainly consider the issue of gap 
anisotropy and competition between different symmetry 
channels in the 2D extended Hubbard model, charac- 
terized by a realistic band dispersion, including nearest 
(N) and next-nearest (NN) neighbors hopping within the 
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CuO planes, and a small-range in-plane potential, allow- 
ing for in-site, N and NN neighbors interaction, in the 
presence of pair tunneling between adjacent layers. 

The issue of the competition among symmetries in the 
gap function arising from the superconducting instability 
of an extended Hubbard model at a given bapd idling 
has previously been considered in the Iiteraturejl3il3 and 
has been recently addressed with renewed attention from 
both the theoretical and experimental points of view, in 
connection with the Cooper pair instability problem in 
lattice fermion systems,EII and with the ^sue of material 
specific phenomenology in the cuprates,Q respectively. 

The ILPT mechanism of high-Tc superconductivity has 
been proposed as a possible natural explanation for the 
observed high values of Tc in layered cuprates, as well as 
a number of other more difficult Jiut related aspects of 
their complex phenomenology.E3l23 On the other hand, 
neither the microscopic origin of the in-plane pairing nor 
the nature of the pairs has to be specified. Several un- 
conventional properties of these materials, due to strong 
correlations already in the normal state, support the idea 
of a breakdown of Fermi liquid theory. In particular, the 
absence of a Drude peak in the low frequency nornial 
state c-axis-aptical conductivity, as observed in YBCOEll 
and LSCOJ23 would rule out metallic transport along the 
c-axis in the cuprates. As a consequence, it has been sug- 
gested that coherent single-particle interlayer tunneling 
is supp f case d- due to the Anderson orthogonality catas- 
tropheJijO'Eil whereas coherent pair tunneling in the su- 
perconducting phase is not restricted. 

Among the mechanisms which would Jirevent single- 
particle tunneling, spin-charge separationE£l has been pro- 
posed. The tunneling process of one fermion would in fact 
require hopping of both spin and charge degrees of free- 
dom, whereas a singlet object, like a Cooper pair, would 
carry charge 2e but no spin. 

Therefore, coherent pair tunneling does not suffer from 
such restrictions, and enters the total Hamiltonian as a 
second order effect in the single particle hopping ma- 
trix element, t^(k), whose dependence on the in-plane 
wave-vector k (see Ref. 19 and Eq. (0) below) has re- 
centlji-been confirmed by detailed band structure calcula- 
tions.CJ The main aspect of the ILPT mechanism is that 
Josephson tunneling of Cooper pairs between adjacent 
CuO layers dramatically amplifies the superconducting 
pairing within each layer, by accessing the normal-state 
frustrated c-axis kinetic energy. 

The addition of such a term to the total Hamiltonian 
does not only greatly enhance Tc, but has also been able 
to describe the observed absence of thje-jHebel-Slichter 
coherence peak in NMR relaxation rate,E2l as well as the 
recent neutron scattering experiments in optirnallv-doped 
YBCOJHa It was also recognized some time agoESTEil that, 
in the same way as the ILPT mechanism very efficiently 
boosts the magnitude of Tc arising from the incipient 
pairing within the planes, essentially due to its near k- 
space diagonality, the amplitude as well as the maximum 
value of the gap function are also dominated by the effec- 



tive coupling induced by the ILPT mechanism. Its actual 
transformation properties under the symmetry operations 
of C^v are however governed exclusively by the intralayer 
contribution to the pairing kernel. 

In this paper, we shall make the latter statement more 
quantitative, showing how the interlayer coupling deter- 
mines the detailed k-dependence of the gap, and actu- 
ally tends to stabilize one symmetry channel compared 
to other possible ones, as the chemical potential is varied 
within the band. 

This paper is organized as follows. In Section || we 
introduce our model and review the basic formalism em- 
ployed to derive the gap equations. In Section |l| we 
discuss the nontrivial numerical problems arising from 
the solution of the latter equations, due to the presence 
of the k-diagonal effective interlayer interaction. A full 
discussion of the gap symmetry structure, its inherent 
anisotropy, its maximum values and locations thereof will 
be included. Calculations of the superconducting density 
of states (DOS) reveal remarkable structures, due to the 
gap anisotropy, which are believed to be relevant to the 
observed anomaloup-phcnomenology in tunneling junc- 
tion characteristics^ In Section IV we shall address the 



issue of determining the critical temperature, as well as 
the temperature at which symmetry mixing occurs, as a 
function of the chemical potential. At exactly T = Tc 
the full k-dependence of the gap function will be derived 
analytically, together with the critical exponents of the 
OP. The resulting expression for the gap function in a 
closed form will serve as an evidence for the non-trivial 
anisotropic character and for the symmetry properties of 
the OP already at the critical point. In Section ^ we 
shall consider various thermodynamical quantities in the 
superconducting phase. Particular attention will be de- 
voted to the normalized jump in the specific heat at the 
critical point, which, at variance with the BCS conven- 
tional result, turns out to be a nonuniversal number, due 
to symmetry competition and to the ILPT mechanism. 
In Section VI we summarize our results and present our 
conclusions. 



II. THE MODEL 

A. Hamiltonian 

The model Hamiltonian we are going to consider in the 
following describes tightly bound interacting fcrmions in 
a bilayer complex: 



H 



k'-kcr'-ko- 



E ^kk'CkT^iki'^k'i'^'T' 



where cjj^ (Ckg-) creates (destroys) a fermion on the layer 
i (i = 1,2), with spin projection a along a specified di- 
rection, wave-vector k belonging to the first Brillouin 
zone (IBZ) of a 2D square lattice, and band dispersion 
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fx, measured relative to the chemical potential 



/i. The second term in Eq. (|lD describes an effective pair 
interaction, already restricted to the spin singlet channel 
only, with 



-^C/kk'<5y 



T,/(k)^kk'(l-4 



(2) 



where N is the number of sites in the square lattice, C/kk' 
measures the coupling interaction within each plane, and 
T/(k) is the tunneling matrix element between adjacent 
layers, motivated by Chakravarty et a/.Ej Equation (|^) 
shows, in particular, how the tunneling mechanism can 
be equivalently described by an interlayer effective inter- 
action term, although local in k-space. 

The main feature of this model is unusual. Although it 
can be cast in the form of a standard BCS-like effective 
Hamiltonian, the second term in the pairing potential 
arises from frustrated kinetic energy along the c-axis of 
the cuprates, unaccessed in the normal state of the high- 
Tc cuprates. However, it is lowered on going into the 
superconducting state. This is a situation which has no 
counterpart in conventional Fermi-liquid based supercon- 
ductors, where the kinetic energy is enhanced upon going 
into the superconducting state, while being overcompen- 
sated by a reduction in potential energy. In the above 
model superconductivity arises via a diametrically oppo- 
site mechanism: Instead of having the gain in potential 
energy overcome the loss of kinetic energy, it is the gain 
in kinetic energy that is the driving mechanism. There 
is ample experimental evidence that the kinetic energy is 
lowered in the superconducting state of the cuprates. Al- 
though a confusing point has been that extracted values 
of the c-axis penetration length has been consistent with, 
estimates of Anderson based on gain in kinetic energy,Ej 
they have also been consistent with c-axis conductivity 
sum rule arguments ignoring the gain in kinetic energy. 
This is traceable to subtleties in applications of c-axis 
conductivity sum rules in unconventional metals, and a 
nice discussion clearing up this CHicial point has recently 
been presented by Chakravartya _ 

Comparison jjpf band structure calculationa^J with 
ARPES resultda for various cuprates suggest that the 
main hybridized single particle band crossing the Fermi 
level can be correctly described in the case of perfectly 
isotropic crystal symmetry by the tight-binding disper- 
sion relation (a being the lattice step) 



Ek ~ — 2t[cos 



s{kxa) + cos(fcj,a)] -I- At' cos{kxa) cos(fcj,a), 

(3) 



where it has been recognizedE3 that at least nearest- 
neighbors (t > 0) as well as next-nearest neighbors 
(t' > 0) hoppings have to be retained, in order to repro- 
duce the most relevant properties common to the mainly 
2D band structure of the majority of the cuprate com- 
pounds. First and foremost, we have in mind the shape of 
the Fermi surface, but also such features as the Van Hove 
singularity in the density of states (DOS), 



(4) 



at /xvH = — 4^', shifted towards the band bottom with 
respect to the mid-band. We hasten to add that we are 
not in any way implying that the Van Hove singulari- 
ties in the single-particle density of states are important 
features in explaii ii j ng || t] |ie large critical temperatures in 
these compounds. c3E30 General, and it seems to us very 
robust arguments for why the Van Ifeve scenario is not 
viable, has been given by Anderson.c2l These considera- 
tions restrict t' /t ^ 0.5, and imply a flat minimum at the 
r point, which gives rise to a pronounced, though finite, 
peak in the DOS at the band bottom. This band has 
a sm^e-particle DOS which can be cast in closed form 



n{s) 



1 



1 



2ttH VI - bCu 



=K 



^ l-[(b + ^)/2p 



(5) 



for 1(a) + b)/2\ < 1, and zero elsewhere. In Eq. (||) we 
have defined b — —2t'/t, a) — e/{2t), and K(a) is a 
compldje elliptic integral of the first kind, with modu- 
lus ao The DOS has a logarithmic singularity n{e) = 
{2^H)-^Vl~^ [log(8/(p - b\)) + log(yT^)] at e = 
2bt, a finite cusp at the lower band-edge n[e = —2t{2 + 
h)] = [47rt(l + 6)]-\ while the DOS at the upper band- 
edge is given by n[e = 2^(2 - &)] = [47rt(l - &)]-^ These 
features of the DOS are important in stabilizing various 
synmietry channels of the OP as the doping level is var- 
ied. E-3 A value of the nearest neighbors hopping param- 
eter ranging around t = 0.25 eV satisfactorily models 
the band structure and the shape of the fijPijii surface 
of the majority of the cuprate compounds .lJo It is not 
among the main aims of this work to specify_the micro- 
scopic origin of the in-plane potential [/kk'-c3 However, 
any potential with the symmetry of the underlying lat- 
tice may be expanded as a bilinear combination of basis 
functions for the irreducible representations of the cryjS-. 
tal point group, which is C/^^ for the 2D square lattice.c2l 
Assuming a finite-ranged potential, a finite subset of all 
the basis-functions (an infinite orthonormal set) will suf- 
fice. Retaining therefore only on-site, nearest and next- 
nearest neighbors in-plane interactions, and projecting 
out interaction terms in the spin triplet channel, one ob- 
tains the following expression for t/kk' , which is separable 
in k-space: 



C^kk' = E A,,5^(k)5^(k'), 



(6) 



))=0 



where (7o(k) = 1, .gi(k) — i[cos(fcj,a) -I- cos(fcj^a)], 
g2(k) = cos(fcj,a) cos(fcya), g3(k) = i[cos(A:j;a) - 
cos(fcya)], g4(k) = sm(kxa) suY{kya) , and [rj — 
0,1,... 4) are phenomenological effective coupling con- 
stants. One immediately recognizes 50 (k), 51 (k), (72 (k) 
to display (extended) s-wave symmetry, whereas 53 (k) 
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and f/4(k) display d-wave symmetry. In the following, 
we shall assume repulsive on-site and attractive intersite 
couphng parameters (Aq > and Ai, A3 < 0), choosing 
their actual values in order to reproduce the correct order 
of magnitude for the critical temperature and gap maxi- 
mum at T = for the cuprates. Throughout this paper, 
we keep A2 = A4 = 0. 

Monte Carlo simulations support the idea that short- 
range antiferromagnetic fluctuations may produce an at- 
tractive intersite interaction (see Ref. |4j for a review). In 
our work, however, such an interaction is taken as phe- 
nomenological, in the sense that an intersite attraction 
is at least required within an extended Hubbard model 
if one expects a d-wave contribution to the OP from the 
lowest lattice harmonics. Remarkably, a perfectly tetrag- 
onal lattice requires Ai = A3. Therefore, if one looks for 
d-wave coupling, i.e. a contribution from 53 (k) to C/kk', 
then one should be also prepared to competition with ex- 
tended s-wave contributions, coming at least from 52 (k). 

Finally, we assume the local dependence of the inter- 
layer pair tunneling matrix element as Tj(k) = t\(}i)/t, 
i.e. a second-order perturbation in the hopping matrix 
element t±(k) orthogonal to the CUiO layers. Recent de- 
tailed band structure calculationaHj formally confirm the 
originjal choice of functional form made by Chakravarty 
et alM 



ti(k) = ^[cos(fca;a) — cosiky 



(7) 



which was arrived at by inspection of ARPES data com- 
bined with analyticity arguments. 

In particular, k-diagonality expresses conservation of 
the momentum component k|| parallel to the Cu02 plane 
during the hopping process. 

We shall see in the numerical cases below that a fine 
tuning of i_L in the range 0.1 — 0.15 eV is the main ingredi- 
ent to reproduce the observed critical temperatures and 
zero-temperature gap maxima in different compounds. 
Such a range is however consistent with band structure 
calculations of tj^.E^ 



B. Mean- field treatment 

A straightforward mean-field (MF) treatment of the 
totaLHamiltonian Eq. (|^) yields the approximate expres- 
sions 

^MF = E ^k^lcL +Y.^K€Ai + H.c], (8) 



ki 



where the auxiliary complex scalar field (i.e., the gap 
function) 



jk' 

= ^Ef^kk'&L,-r,7(k)&i:(i-<5. 



has been introduced. The gap function for the ith layer is 
thus seen to depend on the pair amplitude = (c!_k|Cjj.|) 
in the same layer, through the intralayer potential ?7kk', 
and on the pair amplitude in the adjacent layer, 6^, 
through the interlayer tunneling amplitude rj(k), which 
acts as an effective potential, local in k-space. 

Equation (|^) explicitly shows that, in general, the in- 
terlayer tunneling mechanism endows the gap function 
with a nontrivial, nonlocal structure in the direction or- 
thogonal to the CuO layers. Such a dependence is of 
course relevant in the more general case of multi-layered 
compounds, and its consequences on T,. have been stud^. 
ied, to some extent, by one of the present authors.ca 
A generalization of the methods of the present work to 
multi-layered systems, below Tc is straightforward, and is 
expected to unveil further features in the gap anisotropy, 
due to the coupling of the gap functions in adjacent lay- 
ers. 

In the case of a simple bilayer (i = 1,2), the simpli- 
fying hypothesis that the pair amplitude as well as 
the in-plane single-particle band dispersion and gap 
function Aj^ do not depend on the layer index i allows us 
to decouple the MF Hamiltonian Eq. (||) into, a sum of 
independent Hamiltonians within each layer:E3 

Hmf - E lE^kcLlcL + E[AkcL\cli + H.c] I . 

J I kcr k J 

(10) 

Standard diagonalization techniques in each layer then 
yield for Ak the BCS-like gap equation at the finite tem- 
perature T: 



Ak = -E^kk'Xk'Ak', 



(11) 



where Xk = (2i?k) ^ tanh(/5£^k/2) denotes the pair sus- 
ceptibility, (3 = {keT)-^, and Ey, = -f |Akp the 
upper band of the quasiparticle gapped spectrum. In 
Eq. (^ the pairing potential C/kk' = ;^C^kk' - T'j(k)(5kk' 
includes the finite-range intralayer as well as the local 
interlayer effective interactions. More explicitly. Equa- 
tion dn] ) reads 

1 1 



-■^ l_T,(k)xkiv|?^'''"^'''^'^'- 

Making use of Eq. (||) for the intralayer potential allows 
us to express the gap function as 



with 



(9) 



= - A^-^Ef^('*')^k'Ak'. 



(13) 



(14) 
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At a generic temperature T, Equation (|T^) does not 
yield immediately the explicit k-dependence of Ak as it 
would in the limit of no interlayer tunneling {Tj — > 0). 
This is due to the unusual prefactor [1 — T/(k)xk]"^, 
which includes |Ak| self-consistently via the pair suscep- 
tibility Xk- However, this prefactor manifestly displays 
s-wave symmetry, since | Ak| enters the pair susceptibility 
Xk via the quasiparticle dispersion E^, which is an eigen- 
value of Hmf, and Tj(k) has s-wave symmetry by itself 
[Eq. (^]. Therefore, the complex parameters A^, which 
weigh the basis functions 5,)(k) in Eq. (^3|), measure the 
contributions from the different symmetry channels to 
the full gap function, Ak, at a given temperature T. We 
emphasize that such parameters are not order parame- 
ters in themselves. Only Ak as a whole serves as an OP 
for the superconductive instability, whose onset temper- 
ature Tc is well defined and unique (see Sec. IV below). 
However, a vanishing value of some of the A^ is a sig- 
nal for the absence of the symmetry contribution which 
they represent to the full gap function. Besides, the set 
of parameters {A^} is not unique. Coefficients of any 
other complete orthonormal set of functions would also 
suffice as a basis for expanding the gap function. The 
above choice of basis functions is convenient, since the 
expansion of the in-plane pairing kernel used in this pa- 
per as bilinear combination of basis functions truncates 
after just five terms (t; = 0, . . .4). However, irrespective 
of the choice of basis functions made, when all contri- 
butions are summed up with appropriate weight-factors, 
the result is unique. Moreover, the use of several pa- 
rameters A^ does not attribute to Ak the structure of 
a multi-component (i.e., vectorial) OP. We reserve the 
use of multi-component OP to situations ejicountered in 



systems such as He and possibly UPta 



.00 



The OP of 



high-Tc cuprates is much simpler, with only an amplitude 
and a phase. On occasions, the set {A^} is referred to, 
incorrectly, as the components of a multi-component OP 
of the cuprates. 

We finally remark that the self-consistent expres- 
sion (|l^) endows Ak with an inherently anisotropic k- 
dependcnce, which is modulated by the k-harmonics 
f;,,(k), displaying explicit symmetry. This is best seen 
in the limit T Tc — 0, where the gap is vanishingly 
small, so that Xk in the right-hand side of Eq. ( p^ ) 
can be approximated by its normal state value, viz. 
x'k = (2$k)~^ tanh(/3cCk/2). Only in such a limit, Equa- 
tion (|l^ already yields the explicit expression for the 
k-dependence of the gap function, as a product of the 
anisotropic prefactor [1 — Tj(k)x^'^]~^ and a superpo- 
sition of k-harmonics, weighted with vanishingly small 
coeffici ents A „. We shall precise the latter statement 



in Sec. |IVB| below, where the full k-dependence of an 
incipiently opening gap function at T — Tc will be de- 
rived analytically. It will be shown in Sec. IV that at 
T ~ Tc the presence of a given symmetry contribution to 
the (just opening) gap function generally excludes mix- 
ing with other symmetries. The latter is possible at lower 
temperatures, due to the highly nonlinear structure of 



the gap equations below the critical temperature, at least 
within a given range of the band filling. The mutual ex- 
clusion of orthogonal symmetries in the gap function at 
T = Tc is a well-kn own, result in the case of nonlocal, 
separable potentials .Bca We therefore recover this result 
also with k-diagonal contributions to the potential, such 
as the interlayer pair-tunneling effective interaction. 



III. 



GAP FUNCTION ANISOTROPY AND 
SYMMETRY 

A. The auxiliary gap parameters A,, 



Substituting A^ from Eq. (|]J) into Eq. (Il3|) yields 



with 



ri' 



^'Inv' = ^Xl^k5r7(k)gr,'(k), 



(15) 



(16) 



where Xk ~ Xk/[1 — r,/(k)xk] clearlsi-.acquires the role 
of a 'renormalized' pair susceptibilityJlj These equations 
are in general coupled transcendental equations for Ak, 
and thus define a highly nonlinear problem. Only at 
T = Tc will the situation simplify considerably, as will 
be discussed below. 

However, once self-consistency has been achieved. 
Equations ([l^) are formally linear and homogeneous in 
the phases of the complex parameters A,,, which are re- 
sponsible for the overall complex phase of Ak, as shown 
by Eq. (p^). Due to symmetry considerations, as re- 
marked in Sec. IV A below. Equation ( [Tsl ) reduces to two 
formally independent sets of equations, with real coef- 
ficients, one for each group of parameters belonging to 
either s- or d-wave symmetry. This means that the com- 
plex parameters A^ belonging to the same symmetry are 
all defined up to a same common phase factor. One 
can therefore speak of a relative phase between s- and 
d-wave contributions. In particular, it follows that there 
cannot be anisotropics in k-space of the phase of the or- 
der parameter, other than the (trivial) one arising from 
eventual relative phase differences between two different 
symmetry contributions. (This justifies the widely used 
terminology of s + id symmetry, for example.) 

Due to the presence of the unusual prefactor [1 — 
T7(k)xk]~^ in Eq. (13), which itself must be determined 
self-consistently by finding Ak, ordinary numerical pro- 
cedures used to solve BCS-likC|gap,cauations in the pres- 
ence of separable potentialaiHo'EJO are not applicable 
to the present case. Therefore, remarkably, the gap pa- 
rameters A^ are not enough to define the gap function 
completely: They yield information only about its over- 
all symmetry, on the degree of admixture of the various 
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symmetry channels in the gap function, and on their rel- 
ative phase. The solution for Ak, therefore, has to be 
obtained iteratively for each wave- vector k of interest. 
The iterative numerical procedure employed to solve the 
gap equations is briefly outlined in App. [A[ 

We can now proceed with the solution of the gap equa- 
tion for each given value of the chemical potential fj, and 
temperature T. We first keep /j- at a fixed value. By 
slowly decreasing the temperature from a relatively high 
value, we observe the appearance of a nontrivial solu- 
tion to the gap equations, A^, at a critical temperature 
Tc = Tc{fi), whose value has been made comparable to 
the critical temperatures observed in the cuprates, by a 
suitable tuning of parameters. This onset is signalled by 
a nonvanishing value of some of the parameters A^ , cor- 
responding to a nonzero contribution of one symmetry 
channel (Fig. 0). We shall later show that only one or- 
thogonal channel (restricting ourselves in this work either 
to s- or to d-wave) can contribute to Ak at T = T^, (see 
Sec. 1^ and Ref. Q). Upon further decreasing T below 
Tc, |A^(T)| increases (Fig. 0). Together with |A,,|, we 
plot in Fig. |l| the maximum value of the gap function 
over the IBZ, 



AM(M;r) = max|Ak(Ai;r)|. 

k 



(17) 



One immediately recognizes that Am is considerably 
enhanced with respect to |A^|, which are representative 
of the values it would have had, in the abse nce of ILPT. 

We shall later show analytically (see Sec. IV B| below) 



that A^ and Am behave like ~ {T^-Tf/^ at T = T^-O, 
as it is expected in any mean-field theory for an OP. 
However, Figure ^ shows that the behavior of some A,, 
as functions of T may soon depart from its critical limit 
close to Tc, depending on the value of the chemical po- 
tential. This is to be contrasted with the dependence on 
temperature of Am, which closely resembles the conven- 
tional one in BCS theory. The unconventional tempera- 
ture dependence of the parameters A^ below Tc directly 
stems from their definition, Eq. (p^). A different choice 
of parameters {A,,} would in general lead to a different 
temperature dependence, except their critical behavior 
at Tc- On the contrary, we expect the result obtained for 
Am to be unique, as its value depends more on the ILPT 
amplitude than on the parametrization employed for the 
symmetry character of the OP. 

Depending on the value of /i, other symmetry channels 
may begin contributing to the full order parameter Ak 
as T decreases. This is signalled by a nonzero value of 
the remaining parameters A^ , and by an enhancement of 
the parameters A^ corresponding to the symmetry chan- 
nel already active, as in the numerical example shown in 
Fig. |l|. The critical exponent with which the new A^'s 
open at the criticaL-temperature is again 1/2, as can be 
shown analyticallyiy. The temperature Tm — Tm{f^) at 
which this happens does not correspond to any new in- 
stability: The system is already a superconductor, with 
massive gauge-fluctuations and a finite superconducting 



coherence length. No remarkable feature is to be ob- 
served in Am as a function of T. Its value depends 
more on the anisotropy induced by the interlayer tun- 
neling mechanism than on the intralayer potential. At 
the mean-field level, the OP enhances its overall ampli- 
tude and its anisotropy character, by allowing pairs to 
condense in more symmetry channels. 

Symmetry mixing is made possible by the nonlinear 
character of the gap equations themselves, which be- 
comes increasingly more relevant as the temperature de- 
creases towards T = 0, given our choice of an extended 
in-plane real-space pairing potential. Such a possibilijtjj 
has been already studied in detail by Spathis et al.jlj 
who used a description in terms of a bifurcation OiLjthe 
gap parameters, and by O'Donovan and CarbotteEj in 
the case of an extended Hubbard model without inter- 
layer pair-tunneling. Consistent results have alsO|-been 
obtained by Otnes and one of the present authorsliZi for 
the Cooper problem in presence of an extended intralayer 
Hubbard potential. 

Inclusion of a k-diagonal interlayer pair-tunneling term 
in such a model preserves this feature. A novel effect of 
the interlayer pair-tunneling is that it strongly influences 
the competition between s- and d-wave symmetry chan- 
nels in the OP, enhancing a dominant symmetry channel 
compared with the subdominant other one. The matrix 
element Tj(k) generally reduces the region of symmetry 
mixing in the (/i, T) phase diagram, as will be discussed 
more in detail in Sec. IV. The reason is that when a gap 
amplitude starts to grow at T = Tc, the dominant chan- 
nel will initially suppress pairing in other channels. This 
is generic to any superconductor allowing mixed symme- 
tries to appear in th OP, also conventional ones. Further- 
more, it is important to note that the gap at a certain 
k-point in the BZ depends on the gap at all other k-points 
via the nonlocality of the intralayer part of the pairing 
kernel, even though the inter-layer part is local. The re- 
sult is a strong enhancement of the gap amplitudes in the 
dominant channels, which then to an even stronger de- 
gree will suppress competing channels. The consequences 
of a possible symmetry mixing on observables will be an- 
alyzed below in Sec. ^ We note, however, that due to 
the specific choice of band structure and intra-plane cou- 
pling constants, it is well established that c?-wave pairing 
will dominate in the vicinity of half-filling, while s-wave 
pairing wins out for low filling fractions. Hence, in the 
cuprates, Tj will tend to stabilize d-wave pairing com- 
pared to competing channels, were d-wave pairing to be 
the dominant intralayer channel. 



B. The order parameter |Ak 

Primarily, the interlayer tunneling amplitude Tj(k) in 
Eq. ( p^ ) affects the overall anisotropic structure of the 
gap function, and not its symmetry character. To show 
this, the dependence of A,, on T at a given chemical 
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potential does not suffice alone. Therefore, in Fig. ^ 
we show the overall k-dependence of | Ak| over the whole 
IBZ at T = 0, for a fixed value of fx. We choose to 
plot |Ak| along the family of mutually orthogonal lines 
defined by Ek = const and 7k = const, where 7k is a 
harmonic conjugate of Sk- Such a choice is best suited to 
exhibit and highlight the structure of maxima in the gap 
function along ^k = 0. 

From the numerical analysis, one clearly observes a 
nodal line along the = ky direction for Tm < T < Tc, 
which evolves into a line of local minima as symme- 
tries mix below Tm down to T = 0. Moreover, what 
is more apparent is the presence of rather pronounced 
lines of maxima whose location in the IBZ follow the 
locus of the dispersionless wave-vectors for the normal 
state quasiparticles, i.e. the would-be Fermi line, defined 
by Ck = 0. Absolute maxima (sharp peaks) are located 
at the intersection of the ^k = locus with ky = for 
M < MvH (corresponding to a Fermi line closed around 
the F point), or with k^ — n/a for fj, > //yn (correspond- 
ing to a Fermi line closed around M = (tt/c, tt/o)). Such 
features are of course produced by the enhancing prefac- 
tor [1 — r,/(k)xk]~^ in Eq. (p^, which gives its maximum 
contribtition where r7(k)xk ~ 1, i-e. exactly as quoted 
above. e3 The reason for the sharpness of these features is 
the k-diagonality of the inter-layer pair-tunneling term. 
Similar spikes are difficult to obtain with more conven- 
tional, rie. k-nondiagonal, contributions to the pairing 
kernel ;E3 in such cases we do not get the unusual en- 
hancement factor [1 — Tj(k)xk]~^ in the effective pairing 
susceptibility responsible for the peaks, and anisotropics 
in the pairing kernel tend to be smeared by integrations. 
We suggest that improved energy resolution in ARPES 
is a useful tool to look for sharp features in the gap on 
the Fermi surface, which appears to be a hallmark of the 
ILPT mechanism. 

The maxima distribution and values of Ak along the 
Fermi line is in qualitative and quantitative agreement 
with high-resolution photoemission data available for the 
bilayer Bi2212.l3 It is gratifying to recover such results, 
without making detailed reference to the bilayer band 
structure.E3 It requires invoking the ILTP mechanism, 
where the amplitude Tj(k) depends on k through Eq. (|^) 
(Ref. |l9|) in a-.way which is confirmed by band structure 
calculations .o 

Together with a remarkable k-dependence of the order 
parameter, one observes a different temperature variation 
of Ak depending on the location of k in the IBZ, and par- 
ticularly along the Fermi line ^k = 0, where anisotropy 
is enhanced. This is in qualitative agreement with recent 
ARPES measurements of Ak in-onderdoped Bi2212 at 
different points of the Fermi lineO 



C. Superconducting DOS 

One consequence of such a peculiar anisotropy is e.g. 
^iven by the superconductive density of states at T = 0, 



ns(w) 



where 



Sk 



(18) 

(19a) 
(19b) 



are the usual expressions for the coherence factors in 
BCS-like theories,Ej which hold for an interactiag Fermi 
liquid, in the absence of spectral anomalies.E2l Equa- 
tion ( |l^ ) obviously reduces to n{uj), Eq. (Q), in the hmit 
Ak 0. In Fig. |[ we plot the superconducting DOS 
ns{uj) corresponding to a superconducting spectrum E'k 
with fully anisotropic, prevalently d-wave gap function 
Ak, obtained self-consistently at T = 0, and the analo- 
gous quantity ng{uj), where a pure d-wave gap function 
— A'^gsik) has been used, with A'' = maxk |Ak|. 

In both cases, a gap opens in the SC spectrum at a; = 
(i.e., around the Fermi level). However, the minimum at 

= in ni;{u!) is flatter than in ns{oj), and the fea- 
tures around lu = are quite less pronounced and less 
asymmetric with respect to the Fermi level. Such behav- 
ior in the superconducting DOS is peculiar to-the inter- 
layer tunneling mechanism, and is promisingE3 in order 
to explain the anomalous features observed in tunneling 
junctions experiments with Bi2212.L3 

To complete our picture of the competition of gap sym- 
metries and anisotropy in the ILPT mechanism, we eval- 
uated Ak at T = for chemical potential /i ranging from 
the bottom to the top of the band. In Fig. ^, we plot 
|A^(/i; r = 0)1 against fi. One observes that s-wave sym- 
metry prevails at low band filling, and d-wave symme- 
try at hiahep-fiUing, which is consistent with earlier re- 
sults. olijI^O In a rather narrow region, an OP with 
mixed symmetry occurs. Numerical analysis revealed 
that the ILPT mechanism reduces the extension of the 
latter with respect to the limit Tj — > 0, thus showing that 
a local non-separable contribution to the pairing poten- 
tial frustrates, in general, the coexistence of orthogonal 
symmetries at low temperatures. We argue, therefore, 
that a true, generally non-separable potential, of which 
Eq. (^) is only a truncated expansion over a reduced set 
of basis functions, could even suppress symmetry mixing 
entirely. We have however no formal proof of a such a 
statement, at present. 

From a numerical analysis of the gap maximum at 
T = 0, A?.(/z) = Am{h,T = 0) [Eq. (|^], as a function 
of /i (Fig. Q), we moreover conclude that the ILPT mecha- 
nism yields reasonably large values of thcrgap maximum, 
as observed experimentally in the HTCS,Ej and that the 
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actual values of the intralayer coupling constants con- 
tribute only in a minor way. Furthermore, Fig. ^ shows 
that the largest gaps correspond to prevalently d-wave 
symmetry, and are obtained for fi « /iyn (the exact lo- 
cation depending weakly on A,,), where the enhancement 
due to Tj(k) is highest, once more showing the relevance 
of the 2D character of the single particle dynamics in the 
normal state through their dispersion relation, and the 
importance of the actual value of the next-nearest neigh- 
bors hopping amplitude t' in Eq. (^), which fixes the 
value of at which the Fermi line changes its topology. 



IV. CRITICAL TEMPERATURE 

Among the many experimental facts concerning the 
HTCS phenomenology that the ILPT mechanism is able 
to describe, probably the most apparent is the ease 
with which the high value of the critical temperature 
itself is explained. This is first and foremost due to 
the k-diagonality of the intralayer part of the kernel, 
and has previouslji-been investigated in some detail by 
Chakravarty et al.,i3 when considering the ILPT mecha- 
nism for bilayer compounds such as Bi2212. Of course, it 
is a matter of some importance to investigate the effect 
of inelastic scattering, i.e. k-space broadening, of the in- 
ter layer term, to investigate, how detrimental effect it has 
on Tc ■ Preliminary resultsEa show that Tc is fairly robust 
to a broadening of the interlayer term. 

In this Section, we generalize the results of Ref. |l^ 
to arbitrary doping, conveniently reparametrized by the 
chemical potential ranging within the dispersion band- 
width, extending the analysis to the case of the intralayer 
potential proposed in Eq. (^. The dependence of Tc on 

is a relevant point in itself, since it allows to clarify the 
role of the 2D hole dynamics and that of the incoherent, 
interlayer pair-tunneling mechanism in determining the 
shape and extension of the (/i, T) region allowed for the 
superconductive instability to occur. 

A separate question, in the present context, concerns 
the (/i, T) region allowed to superconductivity character- 
ized by a symmetry order parameter. Due to the struc- 
ture of the gap equation (|T^), such a question involves 
considerable numerical difficulties, in comparisan with 
previous work of some of the present authorsj^ which 
will be dealt with in some detail. 



A. Superconducting instability: pure symmetry 

At T = the mean-field gap function Ak is van- 
ishingly small everywhere in the k-space. Therefore, 
Eq. (nS) linearizes to 



hm M, 



+0 



vv' 



5](V+A,A'C)A,'-0, 

r,' 

where the linearized matrix elements 



(20) 



-y- 

N ^ 1 



^5,(k)v(k) (21) 



do not depend on any more. These matrix elements 
are analogs of the well-known logarithmically divergepl 
integrated pairing susceptibility in the BCS-theory.EHl 
Here, what appears are integrated, effective pairing sus- 
ceptibilities, projected down on various symmetry chan- 
nels. Symmetry dictates that only basis functions having 
the same transformation properties, albeit belonging to 
different irreducible representations of C^y, can yield a 
finite effective pairing susceptibility M°^, . 

The condition for Eq. (|2^) to have a nontrivial solution 
{A^/} is that 



det(5^„, + A^M" ,) = 0. 



(22) 



Due to the s-wave symmetry character of Xk ^^"^ 
definite symmetry character of the basis functions g,, (k) , 
M^^, is block-diagonal. Its elements are nonzero if and 
only if rj and 77' denote symmetry channels belonging to 
the same irreducible representation of the crystal point 
group. Therefore, Eq. (|2^) for Tc at a given fj, factorizes 
into 



T = Tc. 



(23) 



Here, i?^,(M,T,) = det(,5„,,. + KM^^^,) {h = s,d) 

depends only on a subset of the A,, (ry = rjh). Lineariza- 
tion therefore decouples the two symmetries at T — Tc- 
The solution correponding to the largest value of Tc from 
Eq. ( p3| ) corresponds to the true superconducting transi- 
tion temperature. The transformation properties of the 
corresponding eigenvectors determine in which (single) 
symmetry channel the dominant superconducting insta- 
bility occurs. At T = Tc, the other solution correponds to 
a subdominant superconducting instability, and is physi- 
cally irrelevant. Generically, precisely ai T = Tc, we thus 
cannot have an instability into a mixed state, i.e. a su- 
perconducting instability with eigenvectors having com- 
ponents belonging to different irreducible representations 
of C'iy . A mixing of symmetries can only occur below the 
physical Tc, as discussed more in detail in Sec. [H. The 



exception to this statement occurs when /i is fine-tuned 
such that the zeroes of the d- and s-determinants are 
found at the same temperature. The phase space for this 
to occur is however vanishingly small. Such a result is 
a generalization of a known theorem, which applies to 
purely nonlocal separable extended potentials.L3o The 
generalization has been made possible by the definite 
symmetry character (s-wave) of the effective local poten- 
tial induced by the interlayer tunneling amplitude Tj (k) 
in Eq. (H), and is of course extendible to potentials sup- 
porting|_an arbitrary number of symmetry channels in 
the OP.Eij We emphasize that these statements pertain to 
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the square lattice only. In systems with pronounced ab- 
plane orthorhombicity, such as YBCO, a certain amount 
of mixing is expected on quite general grounds, and is 
indeed inevitable. The underlying lattice point group is 
6*21;, and thus an expansion in terms of basis functions 
for Civ will yield several terms. 

The issue of determining Tc = Tdp) and the sym- 
metry channel in which the instability occurs, proves 
therefore to be equivalent to comparing the two solu- 
tions of Eq. (23). It must be noted, however, that at 
variance with the case of no interlayer tunneling, the 
linearized matrix elements M^^, included in the defini- 
tions of {p, Tc) display a divergent behavior at some 
value Tc — T*{^), due to the presence of the denomi- 
nator 1 — r7(k)xk in the-summand of Eq. (|2|). It has 
already been emphasized^ that the self-consistency con- 
dition Eq. ( p^ for a nonvanishing gap Ak below the true 
Tc prevents the occurrence of such a singularity. The 
singularity is due only to the mathematical artifact of 
extending the definitions of the determinants _D"^ Tc) 
to a domain below their zeroes. This is not physically 
meaningful, since the opening of Ak modifies their very 
definitions. 

The occurrence of such an unusual singularity in the 
integrated effective pairing-susceptibility, projected on 
various symmetry channels, has a physical meaning. It 
shows how the action of an interlayer pair-tunneling 
mechanism bounds the critical temperature from below, 
and therefore enhances it. Given an intraplane contribu- 
tion to the pairing-kernel, a lower bound on Tc is set by 
the matrix clement Tj; the lower bound roughly given by 
Tj/4 (see Eq. ( p4| ) below). Beyond this, the actual value 
of Tc is fixed by the intralayer coupling symmetry and 
strength. As already noted, Xu maximum along the 
Fermi line, where limT^T^ lini^fc-+o Xk ~ /3c/4. There- 
fore, the renormalized susceptibility Xk along the Fermi 
line is maximum where Tj(k) is maximum, i.e. at the 
intersection of the Fermi line with the T-X-M path in 
the IBZ (and symmetry related points). Looking for the 
highest temperature T*(^) at which the maximum of Xk 
diverges, one has therefore to distinguish between the two 
possible topologies for the Fermi line arising from Eq. (|^). 
One finds, analytically. 



kBT*{ii) 



64 



\4 

^1 . 



Tj_ ( MT-At 
64 \)iT-2t 



M± < /i < ^iVH, 
/^VH < M < MT, 



(24) 



where , denote the bottom and the top of the band, 
respectively, which generalizes the expression given in 
Ref. At the Van Hove singularity, T*{^) is maxi- 
mum, with fc_BT*(^vH) — Tj/4: — t\/4:t, yielding a lower 
bound keTc < 0.01 eV {Tc < 110 K), which is a repre- 
sentative value for most bilayer cuprates. 

Figure || shows our results for T* and Tc as functions 
of II. The values of the parameters have been chosen as 
quoted in order to yield a critical temperature at opti- 
mal doping whose value is representative of the bilayer 



cuprate superconductors. Superconductivity appears re- 
stricted predominantly to the lower part of the band, 
even though a nonvanishing lower bound T*{p) assures 
a nonzero, albeit decreasing, Tc, as /i increases towards 
the top of the band. In that regime, however, we showed 
numerically that A^pis-aanishingly small (cf. Fig. As 
previously observed s-wave symmetry prevails near 
the bottom of the band, whereas d-wave symmetry wins 
out as fi increases. A robust qualitative argument for 
this was given in Ref. |l^. The critical temperature Tc 
reaches its optimal value near /i = /ivH, the exact loca- 
tion depending on the set of values {A,,} actually chosen 
for the intralayer coupling parameters. 



B. Gap anisotropy at the critical point 

The ILPT mechanism is seen to strongly enhance the 
k-space anisotropy of the gap function also at T = Tc, 
regardless of the symmetry character that the OP takes 
on, which at the critical point is unambiguosly defined 
(no mixing). This is already apparent from Eq. (p^), 
and can be proved by exhibiting the full analytical k- 
dependence of the gap function Ak, at T = Tc — 0. 

For T Tc, one may Taylor expand all quantities of 
interest in powers of /3^|Akp <C 1, safely retaining the 
first nonzero term only. From Eq. {nW, one obtains: 



where 



Ml 



0(/3ek/2) 

[l-Tjik)xi? 



32x3 



(tanh ; 



: sech^ : 



g,(k)v(k)|AkP 
(25) 



(26) 



and a superscript zero denotes that the limit Ak — > 
has been taken. At T = Tc, only one symmetry chan- 
nel is active, therefore Equation (^3|) is satisfied by the 
vanishing of one block determinant, say 13^ (/i, Tc) — 0, 
[h — sovd). Expanding Dx^{fi,T) around /3^|Akp — 
and making use of Eq. (Ea), one finds: 



Dxd^^,T) = DU,,,T)-p^- 



1 



N ^ \1 

k ^ 



mj2) 
Tjik)xl? 



X I J2 A,.9,(k)I^o^,g,,(k) I |Ak| 

, ')')' 



(27) 



+ 0(/3lAk| 

where Wrjr,' denotes the cofactor for the element (5,,,,' + 
XriMriri' in D\^, and a superscript h restricts the sum to 
T] and rj' corresponding to the /i-wave channel only. We 
observe, then, that close to Tc, Equations ( ^ ) factor- 
ize into two separate, independent sets of linear homo- 
geneous equations for the parameters A^, representing 
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either symmetries, respectively. In the proximity of Tc, 
therefore, since (/X, T) = at T = Tc(Ai), only the set 
of equations for corresponding to the incipient /i-wave 
channel admits a nontrivial solution, readily given by: 



|A„| = W°^e, 



(28) 



where fj e {0, 1, 2}, ii h — s, or fj ^ {3, 4}, ii h — d, and 
e is an homogeneity factor, common for all 77's, which 
vanishes as T —> Tc — 0, as specified in the following. In 
deriving Eq. (p8|), we made use of the fact that, in the 
absence of symmetry mixing, such as at the critical point, 
all the belonging to a given symmetry channel share 
the same complex phase factor. Inserting Eq. (|l^) into 
Eq. (p7|), and making use of Eq. (p8|), one finds: 



1 DlM^T) 



(29) 



where 



0(/3Ck/2) 



k ^ 



X I J2 A„5,(k)W^,V5v(k) 

/ \ 2 



X 



(30) 



We finally observe that, by construc- 
tion, limT^T, DUp,T) = 0. Therefore, the expansion 
of Eq. (^9|) around T ^ Tc begins from the linear term 
in (T — Tc), and one straightforwardly obtains: 



a.f (1 



t: 



1/2 



where 



ttQ 

rjrj' 



ly .g.(k)Mk) , 



1 



(31) 

(32a) 
(32b) 



and a superscript c denotes that the limit T ^ Tc has 
been taken. Making use of Eq. ( 31 ) in the expansions for 
|A^| and |Ak|, Eqs. ( p8| ) and (12), respectively, at the 
critical point, one explicitly obtains: 



Tc 



ah- 



T 

T-. 



1/2 



Oc 
h 



lAl 



Tc / T 
""'YV^Tc 



V 



1 - Tj{\,)^i^ 



(33a) 



(33b) 



Equation (33b) analytically yields the k-dependence 
of the gap function at the critical point. In order to 
exhibit more clearly the role of the ILPT amplitude Tj (k) 
in establishing such dependence, one may consider the 
limiting case in which only one basis function (say, 77 = ★) 
contributes to the expansion of Ak. On taking the limit 
jA^/A^J ^ 00, V?7 7^ one recovers the result (see also 
Ref. 



1/2 



.9*(k) 



1-Tj(k)x0^' 



(34) 



where 



1 



(35) 



From Eqs. (^), one also recovers the critical expo- 
nent 1/2 analytically, which is typical for an order pa- 
rameter at the critical poi nt, w ithin a mean-field the- 
ory. Moreover, Equation (33b) clearly shows that no 
symmetry mixing is allowed at T = Tc, by explicitly ex- 
hibiting which basis functions 5,7 (k) contribute to Ak, 
and their weights. The role of the ILPT mechanism is 
furthermore made evid ent b y the presence of the factor 
[l-Tj{\i)x^]-^ in Eq. ( p^ ). This provides the gap func- 
tion Ak with a remarkable anisotropy already at T = Tc, 
thus showing that such an anisotropy is neither due to 
self-consistency (at T = Tc, the values of Ak at differ- 
ent points in the BZ are independent of each other), nor 
to nonlinearity (at T = Tc, the gap equations can be 
linearized). On the contrary, gap anisotropy is robust 
against both self-consistency and nonlinearity, whose rel- 
evance increases as T decreases, as our numerical study 
below Tc has demonstrated. 

From Eq. (33b) one is able to predict a line of rel- 
ative maxima for |Ak| along the ^k = locus already 
at T — Tc- Absolute maxima occur at the intersection 
of the ^k = locus with the T-X-M path, and sym- 
metry related points. The sharpness of the maxima is 
guaranteed by Tj(k), and is therefore distinctive of the 
ILPT mechanism. Away from ^k = 0, the gap function 
is rapidly vanishing over the rest of the IBZ, as an effect 
of the renormalization of the pair susceptibility, induced 
by the ILPT mechanism. Moreover, moving along the 
^k = line in k-space, the gap function |Ak| is seen to 
decrease more than linearly as one approaches = ky, 
where |Ak| attains a minimum value, which is finite and 
very small, in the s-wave case, or zero, in the d-wave case. 
This has to be contrasted with the case of a conventional 
d-wave gap, Ak c>c (73 (k). In such a limit (corresponding 
to Tj —> in our model), Ak would vanish linearly as k 
approaches perpendicularly the nodal line, k^ = ky. A 
flat minimum (node line) along k^ = ky is indeed stronglst 
suggested by ARPES results for Bi2212 single crystals,El 
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and has been earher proposed as-a "smoking gun" for the 
ILPT mechanism by AndersonS 

The sharp anisotropic features of |Ak| are robust 
against non-hnearity, whose relevance increases as T de- 
creases, as shown by Eqs. ( |T^ ) and (|lj). Correspond- 
ingly, the normal state spectrum at T = Tc gets gapped 
where |Ak| is maximum far more significantly than else- 
where. As T decreases, one can think of the the Fermi 
line as remaining practically ungapped along discon- 
nected arcs of ever smaller length. These arcs shrink 
and eventually collapse into a single point along line, as 
T-> 0. 

Recent ARPES experiments ia underdoped Bi2212 sin- 
gle crystals by Norman et aZ£3 are suggestive of such 
a scenario. A progressive 'erosion' of the Fermi line as 
temperature decreases has been related to the opening of 
an unconventional pseudogap, precursor of the supercon- 
ducting gap which opens at Tc. 

V. APPLICATIONS 

In mean-field theory, the superconducting paired state 
is microscopically fully characterized by the gap function 
Ak, which we now have access to over the whole IBZ 
as a function of band filling and temperature. In this 
section, we will discuss how the solution to the gap equa- 
tions can be employed in the calculation of some specific 
thermodynamic properties of the system. A number of 
physical quantities of interest have previously succesfuUy 
been considered with various such solutions, such as for 
instance the. anomalously large gap anisotropy observed 
in Bi2212ji3 the non-conventional featuses in the NMR 
relaxation rate 1/T\ observed in|YBCO,cZl the variation 
of Tc with the number of layers,Ej the unusual faaiures 
in the neutron scattering rates observed in YBC0,c3 and. 
a possible explanation of the spin-gap, or pseudo-gap.E3 
All of the above quoted calculations utilize the special 
features of the gap that arise as a consequence of the in- 
terlayer pair-tunneling mechanism. In particular, the cal- 
culations of the gap anisotropy, the variation of Tc with 
the size of the unit cell, the neutron scattering peak, and 
the spin-gap utilize the unique and sharp k-space fea- 
tures that arise in the solution to the gap equations due 
to the unusual renormalization of the pairing susceptibil- 

ity, Xk^Xk/[l-rj(k)xk]. 

We choose to consider quantities that have the promise 
of being sensitive to the k-space features of the gap, 
which are relatively readily obtained, and which are pos- 
sible to confront straightforwardly with experiments. In 
the following, we shall mainly focus on the specific heat 
anomalies of the model, although work is currently in 
progress concerning the irtplane coherence length and 
the thermal conductivityEj These quantities are either 
sensitive to the presence of the particular Tj-term in the 
Hamiltonian such as specific heat anomalies, or involve 
the derivative of the gap such as the coherence length 



and the thermal conductivity. Moreover we choose, for 
application to the high- Tc compounds, parameters such 
that the critical temperature at optimum doping is given 
by Tc = 90 K. 

Although several of the properties of the supercon- 
ducting state in the high- Tc compounds which in one 
way or another probe the k-space structure of the gap 
are unusual, the thermodynamics seems to be remark- 
ably similar to ordinary superconductors. This is true 
for instance for the entropy of the system. Is a gap aris- 
ing from an unconventional gap equation like the one 
considered in this paper, giving rise to unsual k-space 
features in Ak, consistent with standard thermodynamic 
results otherwise puprmally associated with conventional 
superconductors?L£l Although not shown here, we have 
calculated these quantitites and found that they are re- 
markably similar to those found in any conventional su- 
perconductor .E3 This is basically because quantities like 
entropy involve a k-space integration over smooth func- 
tions of the gap. The detailed k-space features are then 
washed out and the results are to some extent quite in- 
sensitive to these features in Ak. The same also pertains 
to some extent to quantities like the NMR relaxation 
rate 1/Ti, which exhibits features in its T-dependence 
which are reproducible by a gap with a number of dif- 
ferent symmetries. lZH3 This is perhaps not surprising, as 
this quantity involves a double integration over k-space 
vectors. Another matter altogether is the situation where 
we consider quantities involving only one k-space integra- 
tion, and in addition also k.- derivatives of the gap£j 

A. Specific heat 

In this subsection, we will consider the specific heat 
anomalies of the model. The entijopy per particle in the 
superconducting state is given by£2l 

S^fi, T) = -2kB^ 5][/klog /k + (1 - /k) log(l - /k)], 

k 

(36) 

where /k = [1 + exp(/3i?k)]~^ is the Fermi function eval- 
uated with the superconducting spectrum E\^. 

Differentiating S'{fi,T), Eo. (|3|), with respect to T 
one obtains the specific heatEj 

f)Qs 

Q(m,T)=T— (37) 

Whenever i^k, i.e. Ak, contains discontinuities in its 
temperature-derivative as a function of T, the specific 
heat Eq. ( ^ ) displays a finite peak. This is typical of 
the mean-field approximation, as mentioned above. In 
the presence of a competition between several symme- 
try channels, several such discontinuities may occur, at 
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T = Tc and at T = Tm- However, we expect the height of 
the second peak at T = Tm to be exponentiaUy reduced 
with respect to the peak at T = Tr , due to the presence 
of the hyperboUc secant in Eq. (p7|). 

Making use of the gap equations it is possible to de- 
rive a straightforward expression for E\^dE\^/ df3, valid at 
all T < Tc, which turns out to be linear in (9|Akp/9/3 
(cf. also App. ^). Such a quantity is numerically acces- 
sible, in principle, from the solution to the gap equations. 
Therefore, Equation (^) directly yields the temperature 
dependence of Cy also below T^. However, such depen- 
dence turns out to be conventional, and will not be shown 



here (see also Ref. 17). 

At exactly T — Tc, the knowledge of the k-dependence 
of Ak in a closed form allows instead to study analytically 
the jump in the specific heat, normalized with respect to 
the specific heat in the normal state, Cy, i.e. in the 
absence of the gap, at the same temperature: 



CI}- C{>(m,T,) 



(38) 



Making use of Eq. (33b) corresponding to the opening of 
a generic h-wave symmetry gap {h — s,d), one readily 
obtains: 



IE 

6C?r = — fcsa? — V — sech^ I -/J^^k 

C^- = 2kB^ ^ (^IPciu^ sech^ 

We explicitly observe that only at T = Tc one is able 
to include in Eq. (^ for Cy the analytical expressions 
for Ak and its T-dcrivatives: Numerics are only needed 
in performing the integrations over the IBZ where re- 
quired.llj Employin g the value of Tc = T'c(m) numerically 
obtained as in Sec. IV A (Fig. ||), we are eventually able 
to evaluate the normalized jump SCy /Cy'^ in the specific 
heat at r = Tc, as a function of the chemical potential /i. 
We display our results in Fig. ^, and compare them with 
the conventional result SC^/C^" = 12/[7C(3)] ~ 1.42613, 
derived within the BCS theory for an s-wave, uniform gap 
functionEaa 

We find a remarkable agreement with the BCS limit 
over an extended plateau, corresponding to the s-wave 
region in /i. On the contrary, a considerably lower value 
is obtained, on the average, in the c?-wave region, includ- 
ing optimal doping. On the overall, we are thus able to 
predict a nonuniversal ratio 5Cy /Cy'^ , to be contrasted 
with the universai£CS value, valid for Fermi-liquid based 
superconductors.^ This is due to the widely anisotropic 
k-dependence of the gap function (also close to the criti- 
cal point), which is mainly traceable to the renormaliza- 
tion of the pairing susceptibility, and is thus a manifesta- 
tion of the special nature of the interlayer pair tunneling 
mechanism. 



One slightly unusual feature is the possible appearance 
of a second peak in the specific heat at low temperatures. 
Such a feature is not found in a superconductor with an 
order parameter transforming exclusively as a single ba- 
sis function for an irreducible representation of the crys- 
tal point group C/^^. The result we find originates from 
the fact that at low temperatures, new symmetry chan- 
nels couple in to the superconducting order, as shown in 
Fig. This leads to a cusp in the specific heat, hut not 
to any new diverging length in the problem. This second 
anomaly in the specific heat therefore does not repre- 
sent a new superconducting phase transition, but merely 
condensation of Cooper pairs into additional symmetry 
channels. The true order parameter of the problem, Ak, 
becomes finite once and for all at T = Tc, and this point 
represents the only zero-field phase transition. This ap- 
pears-to be a widely misunderstood point in the litera- 
ture.O As mentioned previously, the parameters A,, do 
not represent order parameters for this problem. Note 
that the second anomaly in the specific heat, at low tem- 
peratures, is expected to be well captured by mean-field 
theory. It is located outside the critical region of the nor- 
mal metal-superconductor transition, while this is not the 
case for the first anomaly in the specific heat located at 
T — Tc- Therefore, our results for the main anomaly in 
the specific heat, the prominent step- discontinuity at the 
critical point, should be replaced by a near- logarithmic 
m-jc^gularity characteristic of the ?>DXY-vaoAe\L3 This re- 
Tiects the fact that for optimally doped and underdoped 
systems, phase-fluctuations in the problem appears to be 
|Stoong, such that the true superconducting transition oc- 
curs well below the mean- field transition. 

Our main conclusion of this subsection is that in 
slightly overdoped compounds the main normalized spe- 
cific heat anomaly will be mean-field like, but nonuni- 
versal due to the appearance of a renormalized pairing 
susceptibility Xk/[1 — 2^j(k)xk], in contrast to the stan- 
dard BCS result. 



VI. SUMMARY AND CONCLUDING REMARKS 

We have addressed the issue of the mixing of symme- 
try channels in the superconducting order parameter for 
a bilayer superconductor io-the presence of an interlayer 
pair-tunneling mechanism as a possible framework for 
understanding numerous unconventional features exhib- 
ited by the HTCS compounds. Incipient superconductiv- 
ity has been generated within each individual Cu02-layer 
through a Hubbard-like in-plane potential, including pri- 
marily an on-site repulsion and nearest-neighbor inter- 
action, which has been strongly enhanced through the 
inclusion of the interlayer tunneling amplitude T'j(k), as 
suggested by ARPES as well as by detailed band struc- 
ture calculations. 

A mean-field treatment in the bilayer case allowed a 
computation of the k-dependence of the in-plane order 
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parameter Ak- A suitable numerical procedure has been 
devised in order to solve the inherently nonlinear gap 
equations. It has been possible to study the evolution of 
the symmetry character of the gap function versus tem- 
perature and chemical potential, and unveil a competi- 
tion between s-wave and c?-wave character in At- In this 
description, Ak is a single complex scalar order param- 
eter. No multi-component OP has to be claimed for, 
which would imply the existence of 'more' condensates 
with different 'features,' as elsewhere reported in the lit- 
erature.u In particular, no transition of the normal-to- 
superconductor kind is expected when symmetries are 
allowed to mix: The system is already a superconduc- 
tor with an open energy gap, whose structure in k-space 
only evolves, thus allowing pairs to condense into more 
symmetry channels. 

Moreover, a surprisingly anisotropic (s-wave) pattern 
appeared to be modified by the underlying symmetry 
character. This is evidenced by a strongly pronounced 
line of maxima along the Fermi line, ^k = 0, which closely 
reflected the anisotropic k-dependence of the interlayer 
tunneling amplitude. Such a structure is both qualita- 
tively and quantitatively in agreement jvith the available 
ARPES gap measurements in Bi2212Q and recent phe- 
nomenological gap calculations starting from the multi- 
band structure of the bilayer compounds.EJ Moreover, 
it is essentially embedded in the k-dependence of Tj(k), 
whereas suitable tuning of the intralayer coupling param- 
eters can produce an s-wave contribution which shifts the 
nodes of the gap function slightly away from the T~M di- 
rection {kx = ky), as reported for ARPIj^ experiments 
in bilayer Bi2212 at a given hole content. Q 

The gap obtained within the interlayer pair-tunncling 
mechanism appears to us to be quite promising in ex- 
plaining a number of unusual properties of the supercon- 
ducting state, such as for instance the anpmaJous tunnel- 
ing response observed in HTCS j unctions E!1'E3 Within the 
present approach, such unusal properties are associated 
with sharp k-space features of the gap due to the presence 
of the renormalized pair-susceptibility Xk/[1 — T',/(k)xk]- 

The role of the interlayer tunneling mechanism in en- 
hancing the value of the critical temperature for the 
normal-to-superconducting instabilityJ13 as produced by 
a purely 2D correlation, has been discussed and gener- 
alized for a general doping level, qualitatively reproduc- 
ing the univepS|al non-monotonic dependence of Tc on the 
hole content .O 

The issue of the competition in the symmetry charac- 
ter of the gap function has been addressed both numer- 
ically and analytically in the context of the interlayer 
pair-tunneling mechanism. We were able to verify that 
in the presence of an interlayer pair-tunneling matrix ele- 
ment, the gap symmetry is pure and cannot be mixed on 
an underlying square lattice, at the critical point. The 
gap symmetry belongs to one of the irreducible repre- 
sentations of C'iv, and cannot be expressed as a linear 
combination of several basis functions of such irreducible 
representations. The one exception to this is when the 



chemical potential is fine-tuned to a value such that acci- 
dental degeneracies occur. At exactly the critical point, 
moreover, the full k-dependence and the critical exponent 
of the OP can be derived analytically, thus exhibiting its 
unconventional k-space sharp structure and symmetry 
properties. 

At temperatures well below Tc, and for certain filling 
fractions, mixing of symetry channels may occur. We 
studied the location and width of the (At,r) region al- 
lowing a mixed symmetry superconducting ground state 
on varying the coupling parameters and interlayer tun- 
neling amplitude. In particular, we recovered preva- 
lence of s-wave (resp., d-wave) symmetry at low (resp., 
high) band filling. This is due to the fact that the 
symmetry of the gap is determined by the dominant 
intralayer pairing symmetry, or equivalently the dom- 
inant intralayer dimensionless coupling constant. The 
"symmetry-projected" single-particle densities of states 
of this problem are such that s-wave coupling constant 
dominates at low band- fillings, while c?- wave -coupling 
constants always dominate close to half-fillingEZl. Both 
s-wave and d-wave symmetries are enhanced by increas- 
ing the inter-site attraction, whereas s-wave supercon- 
ductivity is disfavored by increasing the on-site repulsion. 
However, the DOS argument given above also shows that 
the symmetry of the superconducting order parameter is 
highly dependent on doping. 

Finally, we outlined how the solution to the gap equa- 
tion can be employed in evaluating several quantities of 
interest. In particular, we have focussed our attention 
to entropy and the specific heat anomalies of the model 
at the critical point. The entropy in the superconduct- 
ing state is found to have a temperature variation very 
similar to any conventional superconductor, mainly due 
to the fact that it is given by a k-space integral over 
smooth functions involving the gap. The specific heat 
is found to have two unusual features. Firstly, for cer- 
tain filling fractions, a mixing of symmetries may occur 
at a low temperature T = Tm, leading to an anomaly 
in the specific heat, not associated with any true phase- 
transition. Secondly, there is an anomaly at the super- 
conducting transition T = Tc, for which our mean-field 
description is argued to give a reasonable description on 
the slightly overdoped side. This anomaly is analagous 
to the well-known step-discontinuity found in BCS, but 
in our case shows a novel feature. The normalized dis- 
continuity turns out to be not a universal number, not 
only due to the different possible symmetries at the crit- 
ical point, but also depending on the value of Tj, and is 
thus a manifestation of the unusual pairing kernel in the 
gap equation. 
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APPENDIX A: NUMERICAL SOLUTION OF 
THE GAP EQUATIONS 

For any fixed value of the chemical potential /i and 
temperature T, as well as of the coupling parameters 
A,, and interlayer tunneling amplitude Tj, the gap pa- 
rameters An are randomly initialized, and the nonlinear 
equation (|3|) is solved for Ak, for each wave- vector k 
belonging to a suitably chosen fine mesh over the irre- 
ducible sector of the first Brillouin zone, {k : < fca, < 
TT, < fcy < fca;}. The values of Ak thus obtained are em- 
ployed to evaluate M^^/ through Eq. ( [l6| ) . Equation ( |l5| ) 
eventually defines the values of A^, to be used at the 
successive steps in the iteration procedure. The iterative 
procedure terminates when self-consistency is achieved to 
within a preset tolerance limit in |Ak|. Special care had 
to be used near the nodes of the gap function. We ver- 
ified the stability of the convergence procedure against 
the initial choice of A^ , and also by varying the number 
of k-points in the mesh employed in the integrations. 

High accuracy and a resonably small computation time 
in the integration understood in Eq. ([T6|), and else- 
where in the present paper, is made possible by using 
an adaptation te>i the 2D case of the analytical tetra- 
hedron method.L^I An adaptive routine is suitable, due 
to the rapid variation of Xk in Eq. (|l^) for k belonging 
to the locus defined by ^k = (i.e., the Fermi surface 
for noninteracting electrons) . We carefully checked these 
routines by comparing the nurneEically evaluated DOS 
with available exact expressions. EZIl3 
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An index m denotes the inclusion of the gap function 
without the new /i-wave contribution, and that the limit 
T — > Tm has been taken afterwards. 

The above equa tions reduce to the analogous ones de- 
rived in Sec. llVBl at T and Ak = 0. The latter ex- 



pressions are analytical but not in closed form, since they 
require the knowledge of P{d\Ai,\^ /dT), at T = 0, 
which can be accessed only numerically. However, one ex- 
pects such quantity to be vanishingly small for Tm Tc, 
like in the conventional case. (It would be exactly zero if 

Tm = 0). 



APPENDIX B: GAP PARAMETERS AT THE 
MIXING TEMPERATURE 



Fo llowin g a procedure analogous to that outlined in 
Sec. IV B , one is able to derive a critical exponent 1/2 
also for the gap components A^ which open at the mix- 
ing temperature T = T^, thus endowing the gap func- 
tion with an additional contribution with a generic h- 
wave symmetry, orthogonal to the one already present. 
Such a result is consistent with the conventional case 
(Tj — > 0), and is illustrated by the numerical example 
shown in Fig. 0. The calculations are more involved, al- 
though straightforward, and will not be shown here in 
detail. They must however take into account that a gap 
is already open at T = Tm- Generalizing the notation 



introduced in Sec. IVB, the final result is 
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FIG. 1. Temperature dependence of the gap parame- 
ters |A^| (top), and of the gap maximum Am (bottom), 
at = —0.4850 eV. Chosen values of the in-plane cou- 
pling parameters and of the interlayer tunneling amplitude 
are {Ao, Ai, Az, A3, A4} = {0.01, -0.2125, 0.0, -0.2125, 0.0} eV 
and t±_ — 0.08 eV, respectively, yielding a critical tempera- 
ture Tc ~ 13.4 K, at which Ak opens with d-wave symmetry, 
and a mixing temperature Tm ~ 0.2 K, where Ak acquires 
an s-wave contribution. The inset in the top figure shows 
that A3 displays the expected critical behavior, with critical 
exponent 1/2, only very close to Tc. 



I A^l [meV] 




FIG. 2. Dependence of |Ak| over the k in the first quarter 
of the IBZ, at /X = -0.4892 eV and T = 0. Same values of the 
parameters as in Fig. ^ Notice the maxima structure along 
the ^k = locus, including peaks at the intersection thereof 
with ky — and kx = and symmetry related points, 

whose height is enhanced as T decreases, due to the prefactor 
[1-Tj(k)xk]-' in Eq. (|l|). 
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FIG. 3. Superconducting DOS ns{uj), corresponding to an 
anisotropic k-dependent gap in the presence of ILPT (contin- 
uous line), and n|(tj), corresponding to a purely d-wave gap, 
without ILPT (dashed line), at /x = -0.47 eV, T = K. Same 
values of the parameters as in Fig. hi. 
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FIG. 4. Dependence of the gap parameters |A,, (T — 0) 
and of the gap maximum at zero temperature Am on the 
chemical potential p. Same values of the parameters as in 
Fig. I 



FIG. 5. Lower bound temperature T*{jj,) (dashed line) and 
critical temperature Tc{y?j (continous line), as functions of 
the chemical potential fj,. Same values of the parameters as 
in Fig. §. 
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FIG. 6. Normalized jump 5Cy/Cy'^ in the specific heat at 
T — Tc within the ILPT mechanism, as a function of the 
chemical potential (continuous line). Like in Fig. |l|, we used 
{Ao,Ai,A2,A3,A4} = {0.01,-0.2125,0.0,-0.2125,0.0} eV 
and t± — 0.08 eV, respectively. The BCS universal limit 
12/[7("(3)] ~ 1.42613 is also shown, for comparison (dashed 
line) . 
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